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I. INTRODUCTION 

The possibility of altered spontaneous emission by 
modifying the photonic environment, known as the Pur- 
cell effect 0, allows one to tailor the optical emission 
properties of quantum dots (QDs). Spontaneous emis- 
sion is caused by fluctuations of the vacuum electromag- 
netic field, so that its change represents a true quan- 
tum effect. The photonic environment can be altered by 
modifying the density of optical modes, to which the QD 
electronic transitions can couple, and/or by modifying 
the amplitude of the vacuum held at the QD location. 
Both changes have been achieved by placing the QDs in 
a resonator structure with size of the order of the light 
wavelength, in which the electromagnetic field is three- 
dimensionally confined 0, 0, & IE El • As a consequence 
the mode spectrum becomes discretized, and the vacuum 
field amplitude can be significantly modified. 

The light emission of QDs in optical cavities has been 
a very active field of solid state research during recent 
years. The altered spontaneous emission dynamics of 
QDs has been demonstrated using different resonator 
types, such as microdisk structures @, patterned cav- 
ity pillars 0, 0, H, El or photonic crystal defects Q . Its 
experimental verification requires time-resolved photolu- 
minescence (PL) measurements. Since QDs are often 
considered as artificial atoms, it became a standard in 
these studies to use the exponential decay known from 
two-level emitters and to carry it over to a QD system 
in order to quantify the emission dynamics and in par- 
ticular the Purcell effect. Then the decay time r of the 
emitter in the presence of the cavity follows from 
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where tq is the decay time in a spatially homogeneous 
medium, which is determined by the Weisskopf-Wigncr 
decay rate @- The second term on the right hand side 
models the emission into leaky modes. The first term 
describes the QD emission at wavelength A e into a cav- 
ity mode at wavelength A c . An emitter at location r is 
subject to an electric field E(r) whose amplitude varies 
between the maximum value |i? ma 2;| in a field antinode 
and zero for a node position. $ is the angle between the 
electric field vector and the dipole moment of the elec- 
tronic transition. 

The Purcell factor Fp gives the enhancement of the 
emission decay rate in the resonator in comparison to 
the homogeneous medium, 
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Here, Q is the quality factor of the resonator, V c is the ef- 
fective mode volume in the cavity with refractive index n, 
and g is the mode degeneracy. The application of Eq. (1) 
requires that the emitter linewidth AA e is much smaller 
than the cavity linewidth AA C . This is well fulfilled for 
QDs at cryogenic temperatures. 

A closer inspection of the literature reveals, however, 
that in many cases a non-exponential decay of the time- 
resolved PL is observed for a wide variety of QD res- 
onator systems 0, , and even for QDs without 
optical cavities Ll3 • While this in itself complicates 
the quantification of the altered spontaneous emission in 
terms of a constant decay time, we additionally report a 
strong dependence of the time-resolved PL decay on the 
excitation intensity, also far below the stimulated emis- 
sion regime. The decay rate continuously increases from 
the weakest possible pumping, for which we can detect a 
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PL signal, up to the laser threshold. Our experimental 
results are obtained by time-resolved PL measurements 
using (In,Ga)As/GaAs QDs in GaAs-based pillar micro- 
cavities. In these resonators, two distributed Bragg re- 
flector (DBR) mirrors provide the optical confinement 
along the vertical direction. The patterning in the form 
of pillars provides an additional efficient mode confine- 
ment in transverse direction due to total internal reflec- 
tion. Using micro-PL, the emission of individual pillars 
can be analyzed after excitation of the QDs with a short 
laser pulse. 

In the past, various effects have been proposed to ex- 
plain the frequently observed non-exponential decay of 
the PL under these conditions: For example, the exper- 
imental data were recorded from a QD ensemble. As 
the emitter locations vary inside the resonator, for each 
position a different decay rate is expected from Eq. (1). 
Hence, the integrated intensity measured in the experi- 
ment does in general not exhibit a mono-exponential de- 
cay. In addition, the ensemble exhibits fluctuations in 
the QD emission energies and dipole matrix elements, 
which might lead to deviations from a single exponential 
decay. In our calculations, these effects were included. 
For the studied situation, they influence only weakly the 
shape of the time-resolved PL signal. Additional experi- 
mental evidence that the different positions in the cavity 
and variations of the QD emission energies are not the 
prime reason for the non-exponential character of the de- 
cay can be obtained from studies of the time-resolved PL 
signal for QDs without microcavity, which is otherwise 
not the subject of this paper. 01 While the microcav- 
ity enhances the non-exponential shape of the decay, also 
without microcavity the effect can be observed. 

As a further alternative, coupling of bright and dark 
exciton states via spin-flip processes has been suggested 
as an origin for the non-exponentiality 0, 0| . We dis- 
cuss in detail below that this mechanism can be ruled out 
for the present experiments as well. Furthermore, none 
of the alternative mechanisms can explain the observed 
strong dependence on the excitation intensity. 

Based on a microscopic theory for QD carriers inter- 
acting with the quantized light field, we analyze the time- 
resolved PL. The underlying set of equations is referred to 
as semiconductor luminescence equations (SLE), which 
have been used in the past for quantum well systems |l7j . 
While the PL of an (ideal) ensemble of two-level emit- 
ters shows an exponential decay with a time-constant 
independent of the excitation conditions, this is found to 
be very different for semiconductor QDs. The recombi- 
nation of an excited electron requires the presence of a 
hole. Scattering and dephasing processes reduce the cor- 
relation between optically generated electron-hole pairs 
and, thus, lead to a distinct departure from the simple 
two-level picture (corresponding to independent excitons 
with fully correlated electron-hole states). Our calcula- 
tions of the electron- hole recombination reveal the intrin- 
sic nature of the non-exponential decay and its excitation 
intensity-dependence. While other (extrinsic) processes 



could be important in other specific experiments, they 
do not show the discussed intensity dependence and are 
limited to special systems or excitation conditions. 

The article is organized as follows. In the next sec- 
tion we discuss the samples under study as well as the 
experimental technique. In Section 3 we present the ex- 
perimental results and give a preliminary analysis. In 
Section 4 the theoretical model is introduced and the re- 
lation to the experiment is discussed in detail. The article 
is concluded by a summary. 



II. SAMPLE AND EXPERIMENT 

The planar microcavity sample was grown by molec- 
ular beam epitaxy on a (lOO)-oriented undoped GaAs 
substrate, with a GaAs buffer layer of 0.4 /jm thickness. 
The GaAs A-cavity layer was sandwiched between two 
DBRs, consisting of 23 and 20 alternating AlAs/GaAs 
films for the bottom and the top mirrors, respectively. 
Each film is made from a 79 nm-thick AlAs and a 67 nm- 
thick GaAs A/4 layer. A single layer of self-assembled 
(In,Ga)As/GaAs QDs serves as optically active medium 
in the center of the resonator, where the vertical electric 
field amplitude has an antinode, to maximize the light- 
matter interaction in the planar cavity case. The nomi- 
nal material composition of the QDs is InAs, but during 
growth intermixing with the GaAs barriers occurs. As 
the precise, position-dependent material composition is 
not available, we use the generic term (In,Ga)As for the 
QD material. The QD surface density is w 3 x 10 10 cm~ 2 . 
Single pillar microcavities with different diameters rang- 
ing from about 1 to 6 fim and spaced 400 fim apart were 
fabricated by electron beam lithography and dry etching 
[l8j . Due to this patterning the field strength is modu- 
lated in the cavity plane. 

The pillar microcavities were mounted on the cold fin- 
ger of a microscopy flow-cryostat allowing for tempera- 
ture variations down to 6 K. In time-integrated photolu- 
minescence spectroscopy, a frequency doubled Nd:YAG 
laser was used for continuous wave optical excitation at 
A = 532 nm. For time-resolved photoluminescence, a 
pulsed Ti:Sapphirc laser with pulse durations of wfOO fs 
was used. The laser beams were focussed onto the sam- 
ple by a microscope objective with a focal length of 1.3 
cm, by which a spot diameter of about 10 /jm could be 
reached. This diameter is larger than that of the largest 
studied cavity, so that one can assume homogeneous ex- 
citation conditions. 

We note, that the laser excitation does not lead to any 
sizeable sample heating effects in the microscopy cryo- 
stat. This has been checked by performing time-resolved 
experiments with the sample inserted in superfluid he- 
lium (T = 2 K) in an optical bath cryostat. Under oth- 
erwise identical excitation conditions, the same behavior 
is observed for the micropillar photoluminescence kinet- 
ics as the one described in Section III. 

In particular, the measured time evolution of the PL 
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depends strongly on the excitation power. In the fol- 
lowing we give the average power density P exc , which is 
connected to the energy per laser pulse J pu ise, focussed 
into a spot with area A, by: P eX c X A = f x J pu ise, where 
/ = 75.6 MHz is the pulse repetition rate. A laser power 
of 1 mW (corresponding to a power density of 1.27 kW 
cm~ 2 , which is a typical value for high excitation in the 
experiment) is achieved for a pulse energy of ~ 0.013 
nJ. With an excitation energy of 1.4 eV per electron- 
hole pair, about 5.8 -10 8 of such pairs could be created. 
Further we assume that within the A-cavity only about 
0.01% of the incident laser power is converted into carri- 
ers, which may be trapped in the wetting layer and relax 
further towards the QD ground state. This is estimated 
from comparing the relative emission intensities of the 
QDs with that of bulk GaAs. Distributing 0.01% of the 
carriers over the excitation area leads to an estimate of 
n e h ~ 7.5 x 10 9 cm -2 for the carrier density. Carrier den- 
sities in this range have also been used in the numerical 
analysis, see Section llVl 

Some effects have not been considered in this simple 
estimate, because they are difficult to quantify. For ex- 
ample, the true carrier density may be reduced due to 
above-stopband reflection. On the other hand, an in- 
crease might occur due to reabsorption of light emitted 
from GaAs, such as from the substrate. Our cavities op- 
erate in the weak-coupling regime, so that effects known 
from strong coupling play no role here. Note that the 
variations in the carrier density due to these additional 
effects are expected to be small. 

Excitation and collection were done through the micro- 
scope. After light collection, the emission was directed 
into a 0.5 m monochromator where the signal could be 
sent either to a charge coupled devices camera for time- 
integrated PL studies (used also for alignment) or to a 
streak camera for time resolved experiments with a res- 
olution of fts 20 ps. 

In Fig.fjthe emission spectra of single pillars with dif- 
ferent diameters arc shown, obtained by excitation with 
the Nd:YAG laser. For better comparison, the intensity 
has been normalized. With decreasing size the energies 
of the optical modes shift to higher frequencies. In addi- 
tion the splitting between the modes increases strongly. 
These observations are in accordance with previous stud- 
ies of the optical mode spectrum in similar patterned 
cavities [U IH HHH 111 H| . We note that the inho- 
mogencously broadened QD emission spectrum has a full 
width at half maximum of about 30 mcV at low excita- 
tion. The emission is centered around 1.38 eV, so that 
the QD ensemble represents a light source which provides 
emission over the energy range in which the fundamental 
modes of the studied optical resonators are located. 

We have studied the PL decay of the QDs across their 
inhomogeneously broadened emission band (by analyzing 
the emission along the plane of an unpatterned resonator 
and by studying a QDs reference sample). An initial de- 
cay time of 600 ± 50 ps is observed at 0.17 kWcm -2 
with no correlation to the emission energy. Therefore 
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FIG. 1: PL emission spectra from single cylindrically shaped 
micropillar structures of varying diameters. The inset gives 
the quality factors of the two energetically lowest photon 
modes as function of the pillar diameter. Solid (open) sym- 
bols give data for the fundamental (first excited) mode. T = 
6 K. 



any cavity size dependence of the carrier lifetime cannot 
be related to systematic variations of the dipole coupling 
with emission energy due to changes of the QD confine- 
ment. 

The inset of Fig. Ogives the pillar diameter dependence 
of the quality factors Q — E/AE of the optical modes. 
Data for the two lowest confined modes are shown: in 
both cases we find a considerable decrease in Q with de- 
creasing pillar size. While for the 6 fjja diameter pillars 
the quality factors are almost 10000, the cavity quality 
for the smallest cavities with 1 /im diameter varies de- 
pending on the resonator from 2000 to below 1000. Also 
for larger cavity diameters, sample-dependent variations 
of the cavity quality are observed. Furthermore, the first 
excited mode has always a smaller Q than the fundamen- 
tal mode. However, one has to be careful in such a com- 
parison, as an increased line width of the corresponding 
emission might arise from a slight mode splitting, which 
has been predicted in Ref. [2f| for the first excited mode. 
Within our experimental accuracy (0.2 meV resolution 
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FIG. 2: (color online) Low excitation time-resolved PL emis- 
sion for micropillars with different diameters. The excitation 
power density was 1.3 kWcm~ 2 . The decay times correspond- 
ing to the single exponential fits shown by the solid lines are: 
400 ps (6/im), 315 ps (4/xm), 200 ps (3/im), 110 ps (2/xm), 
and 80 ps (1/im). For clarity, the traces have been shifted 
vertically. 



of the setup), this mode splitting lies below the mode 
line width due to the finite photon lifetime. Independent 
of the involved photon mode, the decrease of Q with de- 
creasing pillar size arises from reduced confinement of the 
mode for small cavities and also from an increased im- 
portance of surface roughness scattering at the sidewalls. 



III. 



EXPERIMENTAL RESULTS 



Figure [2] shows the decay of the time-resolved PL for 
detection at the energy of the respective fundamental op- 
tical mode of micropillars with different diameters. The 
intensity is plotted on a logarithmic scale. Excitation was 
done with the pulsed Tisapphire laser, with the wave- 
length set to 800 nm, corresponding to creation of car- 
riers in the GaAs barriers, to allow for a variation of 
excitation power and therefore carrier density over wide 
ranges. This wavelength is also above the stop-band of 
the planar resonator. The used low excitation power of 



1.3 kWcm 2 guarantees that the observed PL occurs in 
the spontaneous emission regime. 

The faster decay for decreasing pillar diameter is 
mainly a consequence of the confinement induced en- 
hancement of the vacuum field amplitude which results 
in the Purccll effect. The reduction of mode volume leads 
to an increase of the Purcell factor, as for the discussed 
range of diameters the mode volume decreases faster than 
the Q-factors (see inset of Fig. ^| . For the decay of the 
signal over the first order of magnitude, the deviation 
from an exponential decay is rather weak. Straight lines 
have been added to fit a Ti time to the initial decay. 
However, on a larger scale the decay data clearly reveal 
a non-exponential character. 

Here we note explicitly that the data cannot be de- 
scribed by biexponential or stretched exponential de- 
cay forms. While such forms naturally can match the 
data better than monocxponential decays, as they involve 
more fit parameters, they still result in considerable de- 
viations from the data. Good agreement can generally 
only be reached by multiexponential decays involving a 
large number of parameters without physical meaning. 
Biexponential forms would be appropriate if there were 
two independent decay channels, each with a consider- 
able contribution to the emission. Potential candidates 
for additional decay channels besides the excitonic one 
such as spin-dark excitons, charged excitons, and so on 
will be explicitly ruled out by the arguments given below. 

To obtain more insight into the emission dynamics, we 
have varied the excitation power density P exc . Fig. |3| 
shows the time-resolved emission of a 5 /im pillar for dif- 
ferent P eX c at an excitation wavelength of 800 nm. With 
increasing power, the decay becomes generally faster. 
Even for the lowest P eX c, for which we could record a 
time-resolved signal, no saturated, power-independent 
decay is found. For the highest excitation powers, rise 
and initial decay become more and more symmetric with 
respect to the signal maximum, indicating that the res- 
onator has been pushed into the stimulated emission 
regime. Similar behaviors have been observed for pillars 
with other diameters. 

First, we need to address the question, to what ex- 
tent the observed dependencies regarding cavity size and 
excitation power are influenced by non-radiative decay 
channels, such as traps at the etched cavity sidewalls. 
Increasing importance of such traps with decreasing cav- 
ity diameter might also lead to the lifetime shortening 
observed in Fig. 2. To analyze this effect we have per- 
formed studies at varying temperatures. For T > 50 K 
thermal emission out of the QD confinement becomes im- 
portant, leading to a drop of PL intensity. At lower T, 
however, the integrated intensity is constant and also the 
decay times do not vary with temperature, which indi- 
cates a negligible influence of non-radiative decay. In this 
regime, the emitter linewidth AA e clearly falls below the 
cavity mode linewidth AA C , which is a prerequirement 
for the observation of the Purcell effect. 

The minor importance of non-radiative decay is also 
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FIG. 3: (color online) PL decay curves of a 5 /im cavity 
at different excitation powers. The decay times correspond- 
ing to the single exponential fits shown by the solid lines 
are: 550 ps (0.17 kWcm" 2 ), 475 ps (0.51 kWcm" 2 ), 265 ps 
(1.7 kWcm" 2 ), 120 ps (5.0 kWcm" 2 ), 30 ps (7.6 kWcm" 2 ), 
and 20 ps (9.1 kWcm~ 2 ). T = 6 K. For clarity, the traces 
have been shifted vertically. 



supported by a variation of the experiment, where the 
excitation wavelength is changed to 860 nm, which is 
slightly above the resonator stop-band edge, but below 
the GaAs barrier and into the wetting layer. This wave- 
length shift affects carrier capture and relaxation, but 
should not influence the cavity size dependence of the 
decay rate. However, it would be of importance if non- 
radiative decay were relevant, as for excitation above the 
barrier the photogenerated carriers may diffuse to the 
cavity sidewalls, while for excitation into the narrow wet- 
ting layer (which is subject to carrier localization effects 
at low temperatures) diffusion to the sidewalls is strongly 
hampered. Our data show no change in the decay rate 
shortening with decreasing cavity size for the two excita- 
tion conditions. The entirety of these tests allows us to 
relate the PL decay time to the radiative decay. 

Different reasons for a deviation of the PL emission dy- 
namics from an exponential decay need to be considered. 
For very low excitation power, contributions from exci- 
ton complexes such as biexcitons can be ruled out. How- 



ever, charged exciton complexes may be formed due to 
unintentional background doping, for example. We have 
tested the presence of residual charge carriers on a QD 
reference sample by Faraday rotation measurements |27j 
in a magnetic field normal to the heterostructure growth 
direction. Presence of free carriers would lead to observa- 
tion of pronounced spin quantum beats which last longer 
than the exciton lifetime. As such beats could not be 
observed, we can safely conclude that the vast majority 
of QDs is undoped. 

Further, in the present case the excitation was non- 
resonant allowing for fast spin relaxation in the barriers. 
Therefore not only spin bright excitons are formed, but 
also spin dark excitons, which can recombine radiatively 
only after a spin-flip. Previous investigations have shown 
that, at cryogenic temperatures, spin-flips are strongly 
suppressed for carriers in the QD ground states. For elec- 
trons the only viable mechanism seems to be the hyper- 
fine interaction with the lattice nuclei, while spin-orbit 
interaction has been shown to give spin-flip rates in the 
kHz-range only, corresponding to much longer time scales 
than the ones considered here psj . Spin-orbit interaction 
is also the only viable mechanism for the holes, but at 
low temperatures a two-phonon process is required to 
induce a spin flip [2^ . In agreement with these consider- 
ations, exciton spin- flip times have been reported which 
are much longer than its radiative decay time |3Cj . 

If feeding of the spin-bright exciton reservoir from the 
reservoir of dark excitons were important, two time-scales 
would be relevant: Apart from the bright exciton decay 
during times < 1 ns the dark exciton background would 
decay on time scales of nanoseconds or longer. From the 
study of a QD reference sample we find that no such 
background can be identified at cryogenic temperatures. 
For the time scales of interest it would appear as contri- 
bution to the constant background due to dark counts, 
which does not affect the decay analysis, as this back- 
ground is subtracted. Thus, also the dark exciton states 
cannot be thought of as origin of the non-exponential 
decay at the used low temperatures. However, a strong 
background with decay times in the few ns range could 
be observed by raising the temperatures to a few tens of 
K, where the thermally induced phonon population can 
lead to flip processes through spin-orbit coupling. This 
is confirmed by the background decay time shortening 
strongly with increasing temperature. 



IV. THEORETICAL MODEL 

The aim of this section is to outline a theoretical model 
for the PL dynamics of QDs, which includes population 
effects in the carrier system, the many-body interaction 
between the carriers as well as their interaction with the 
quantized light field within a microcavity. 

QDs are often compared to atomic systems due to 
the appearance of localized states with discrete energies. 
However, QDs usually contain many electronic states and 
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excitations involve many electrons and holes, which are 
influenced by the Coulomb interaction. Additional carri- 
ers in the WL states contribute to screening and dcphas- 
ing which - together with scattering processes among the 
localized carriers - weakens correlations between elec- 
trons and holes. 

This situation differs fundamentally from light-matter 
interaction of atomic two-level systems, which are often 
used for a simplified analysis of QDs. A two- level model 
is applicable if the optical field couples resonantly only to 
two electronic levels and if the excitation involves only a 
single electron. In this case the appearance of an electron 
in the upper state is inescapably linked to the nonexis- 
tence of an electron in the lower state. In the semicon- 
ductor language, electron and hole populations are fully 
correlated. As a consequence, the time-resolved PL of 
two-level systems shows an exponential decay. 

To describe the QD PL, we use the SLE for a system 
consisting of interacting charge carriers and a quantized 
light field, which has previously been applied to quan- 
tum wells 0. The SLE describe the coupled dynamics 
of the electron- and hole-population, f„' , the general- 
ized photon population, (b q b q i), and the photon-assisted 
polarization, (b q h a e a ), in the incoherent regime: 
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Here, e' a , e a and e e a denote the creation- and annihilation- 
operators and the single-particle energy of an electron in 
state (jf a (r). The corresponding quantities for the holes 
are h/ a , h a , e^, and </>Jj;(r). The states for charge carriers 
are either delocalized WL or localized QD states. The 
operator b q (b q ) creates (destroys) a photon in the opti- 
cal mode q, which is characterized by the complex res- 
onance frequency cu q and the transversal mode-pattern 
u q {f). The light-matter coupling is determined by g qa cx 
/ d 3 r(f>" (f)efu q {r)(jy^{r) . The exchange- Coulomb matrix 
elements are denoted by V a p and the single-particle en- 
ergies including Hartree-Fock renormalizations are given 
by ia' h \ The population changes due to scattering are 
treated in a relaxation time approximation. 

Eqs. © and (0} show that the dynamics of the car- 



riers and photons are driven by the photon assisted po- 
larization, which is governed by Eq. JHJ). The first line 
of Eq. 10) describes the free evolution, the second line 
is responsible for excitonic resonances in the spectrum 
and the third line describes stimulated emission or ab- 
sorption. The last line contains the source term due to 
spontaneous emission, which dominates over the stimu- 
lated term for weak excitations. Nevertheless, in high-Q 
microcavities, the rcabsorption of photons can modify the 
results for the time-resolved emission even for weak exci- 
tation. On the Hartree-Fock level (corresponding to un- 
corrected carriers) the source term is given by g qa fafa- 
Correlations due to Coulomb interaction of carriers are 
included in f2™£ and are evaluated on singlet level |3lj . 
We investigated the influence of higher-order contribu- 
tions (doublet level), as they are discussed in js^l for a 
quantum well, and find that they play no important role 
in the presence of the feedback provided by the cavity. 

To gain additional insight into the physics described 
by the SLE, we use - only for the following discussion 
in this paragraph - some simplifications: We disregard 
stimulated emission/absorption and neglect the Coulomb 
interaction. The adiabatic solution of Eq. then yields 
for the population dynamics 



dt 



J a 



(e,h) 



J a. J a 



opt 



(6) 



where r sp is the time-constant for spontaneous emission. 
The HF contribution to the source term f^fa clearly 
leads to a non-exponential decay. Furthermore, the rate 
of decay depends on the carrier density and is higher 
for larger population. For the calculations presented 
below, none of the mentioned simplifications have been 
made: stimulated emission/absorption are included and 
Coulomb effects such as corrections to the HF factor- 
ization of the source term of spontaneous emission are 
considered. 

The discussed experiments can be used to distinguish 
between the two regimes of fully correlated electron-hole 
pairs, leading to an exponential decay of the PL, and 
partially correlated carriers with a non-exponential PL 
decay. As one is typically interested in the luminescence 
decay over a ns time-scale rather than the initial exci- 
tation and relaxation of the system (which takes place 
on a ps time-scale), we focus here on the dynamics that 
occurs after the system has been excited and the coher- 
ent polarization has decayed due to dephasing. In this 
incoherent regime we can use a Fermi-Dirac distribution 
of carriers in the WL and QD states and zero photons in 
the cavity as our initial conditions and evolve the system 
according to Eqs. JSJ-JSJ. 

In order to correctly account for the effects of the size 
distribution of the QDs (different transition energies), 
their spatial distribution inside the cavity, and the dif- 
ferent dipole orientations (different coupling matrix ele- 
ments), it is not sufficient to analyze the equations for a 
single dot with averaged properties. Instead it is neces- 
sary to solve the SLE for an entire ensemble of different 
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FIG. 4: Calculated PL for an ensemble of QDs in a 6 fim 
diameter pillar microcavity with initial carrier densities from 
1 x 10 9 cm~ 2 to 5 x 10 9 cm~ 2 in equidistant steps from top to 
bottom. For better comparison the results are normalized. 



QDs. For our calculation we take only a fraction of the 
total number of QDs with a transition frequency close to 
that of the relevant cavity modes. Therefore, the effec- 
tive density of QDs resonantly interacting with the cavity 
mode is assumed to be 3 x 10 9 cm -2 , distributed in an 
interval of approximately 1.5 meV. 

The individual QDs are modeled with a harmonic con- 
finement potential in the WL plane and a step-like con- 
finement in growth direction. The strength of the har- 
monic confinement is varied for different QDs to account 
for the inhomogencous broadening typically observed in 
this material system. We restrict our analysis to QDs 
with s- and p-shells for electrons and holes. 

The transverse mode-pattern and resonance frequen- 
cies uj r q cs of the optical modes are calculated using a 
three-dimensional transfer-matrix approach (for details 
see [Iff). The corresponding quality- factors Q are ob- 
tained from the experiment. The complex resonance fre- 
quency is then given by cu q — wj cs (l — i/Q)- While it 
is sufficient to include only one resonant mode for the 
smaller pillar, for the larger pillars several modes have 
to be taken into account. The coupling between differ- 
ent modes is neglected, {b\b q >) k 5 qq >(b q b q ). Besides the 
resonant modes, which are characterized by their large Ca- 
valries and pronounced peak structure in a transmission 
spectrum, there exists a background contribution from 
the continuum of leaky modes. In order to include their 
influence, we assume that the background contribution 
consists of a fraction of the continuum of modes of the 
homogenous space. The size of this fraction can be esti- 
mated by counting the plane waves that i) either reach 
the sidewalls of the micropillar in an angle smaller than 
the critical angle of total internal reflection, or ii) have 
a momentum component fc| along the pillar axis that 
lies outside the stopband of the DBR and can therefore 
immediately escape from the cavity. 

In Fig. 0] we show the number of photons in the funda- 




Time [ps] 

FIG. 5: Calculated PL of QDs in a pillar microcavity with 
various diameters for an initial carrier density of 2 x 10 9 cm -2 . 



mental mode leaving the cavity per unit time. Different 
initial carrier densities are used to model the variation 
of the excitation power in the experiment. The non- 
exponential decay of the signal is clearly evident. Fur- 
thermore, the rapidness of the decay strongly depends on 
the initial carrier density, which corresponds to different 
pump intensities in the experiment. This shows that it is 
not meaningful to introduce a decay time that depends 
only on the photonic density of states without including 
the influence of the carrier system. Instead, a thorough 
analysis of time-resolved PL signals has to take both the 
carrier system and the photonic system into account. 

It should be noted that for strong optical fields the gen- 
erated carrier density depends in a nonlinear manner on 
the pulse intensity due to saturation effects. It is not the 
purpose of the paper to quantify these optical nonlin- 
caritics together with the subsequent carrier relaxation 
and to directly connect experimental pump intensities 
and the resulting carrier densities. Instead we focus on 
the physics of the recombination dynamics and empha- 
size the strong dependence of the time-resolved PL decay 
on the carrier density in the system. 

The calculated PL for fixed initial carrier density but 
different diameters of the micropillar cavity is displayed 
in Fig. [S] The smaller pillars show a faster decay in 
connection with a larger Purcell factor, as has been dis- 
cussed in Sec. IIHI The different heights of the curves 
can mainly be attributed to the fact that in larger pillars 
more carriers distributed over more QDs take part in the 
recombination dynamics. 

Frequently it is argued, that the non-exponential de- 
cay observed in PL measurements stems from a superpo- 
sition of many exponential PL signals of various emitters 
with different cavity positions. The role of an inhomo- 
geneous distribution of QDs is analyzed in Fig. [S] The 
solid line represents the calculated decay of the time- 
resolved PL from an ensemble of QDs with various cav- 
ity positions and fluctuations of the transition energies 
and dipole moments (same as in Figs. 4 and 5). For 
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distribution effects. 

The deviations of the measured from the calculated re- 
sults close to time t = and in particular the somewhat 
slower raise of the PL signal observed in the experiment 
can be attributed partially to the fact that the optical 
carrier generation was not modelled and that the experi- 
mental setup has only a finite time resolution. Note that 
we did not adjust the calculations to have a quantitative 
agreement between experiment and theory. Too many 
parameters, such as size and composition of the QDs, 
are unknown in detail. 



01 1 — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

300 600 900 1200 1500 

Time [ps] 

FIG. 6: Calculated PL of an ensemble of QDs with inhomo- 
geneous broadening (solid line) , for identical QDs with maxi- 
mum coupling strength and on-resonance transitions (dashed 
line), and for identical QDs with averaged coupling strength 
(dotted). The pillar diameter is 6 /J,m and the initial carrier 
density is 2 x 10 9 cm~ 2 . 

comparison, the dotted line shows the result for identi- 
cal QDs with averaged values for mode-strength, cavity 
field, transition energy, and dipole coupling. While the 
decay remains non-exponential, the decay rate is strongly 
underestimated. If identical QDs with maximum values 
for mode-strength, cavity field, and dipole coupling as 
well as resonant transition energies are assumed (dashed 
line), the decay rate is slightly overestimated in the ex- 
ample with practically the same shape as for the inho- 
mogeneous QD distribution. This result shows that the 
QDs with efficient coupling to the cavity field dominate 
the emission properties. The non-exponential character 
of the decay is only weakly determined by inhomogeneous 



V. SUMMARY 

A microscopic description of the QD emission, based 
on the "semiconductor luminescence equations" includ- 
ing many-body Coulomb effects, shows the appearance 
of a non-exponential decay that is intimately connected 
with the intensity dependence of the decay from weak 
to strong excitation conditions. The results explain the 
time-resolved PL of QDs in pillar microcavities. Other 
origins such as contributions from spin-dark excitons, 
charged excitons etc. have been ruled out as dominant 
contribution to the observations for the presented exper- 
iments. 
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